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We generalize the Gutzwiller wavefunction for s = | spin chains to construct a family of wavefunc- 
tions for all s > i. Through numerical simulations, we demonstrate that the spin spin correlation 
\l ' functions for all s decay as a power law with logarithmic corrections. This is done by mapping 

>--^ , the model to a classical statistical mechanical model, which has coupled Ising spin chains with long 

range interactions. The power law exponents are those of the Wess Zumino Witten models with 
k = 2s. Thus these simple wavefunctions reproduce the spin correlations of the family of Hamilto- 
r^ , nians obtained by the Algebraic Bethe Ansatz. 

^■ 
^ ! INTRODUCTION 

(N; 

^ Quantum spin chains with power law correlations are of great current interest. Starting with the Bethe chain 1], it 

^~Ti I has been known for a long time that half integer spin systems generically have spin correlation functions that decay 

I ■ as the inverse of the distance between the spins. For integer spin chains, the correlation functions are generically 

^_> I expected to decay exponentially, due to the Haldane conjecture Q- However, it is also known that special models can 

. ■ be constructed for which the spin correlation functions decay as a power law for all (including integer) spin s, with 

jrt I an exponent that depends on s. These models are realized from the Bethe Ansatz in its algebraic form. The work of 

Cj ■ the Leningrad school has provided concrete realizations, e.g. the model of Takhtajan and Babudjan [2j for spin 1. It 

I ' has also provided a general technique for obtaining the Hamiltonian for all values of spin, starting from the Algebraic 

^ . Bethe Ansatz 3- 

5 The coefficients of the biquadratic and other terms need to be specially chosen to obtain these models. Alternatively, 

the Algebraic Bethe Ansatz method automatically generates these Hamiltonians. The general construction of a field 

theory for these models was undertaken by Affleck pj, who showed that these models belong to the Wess Zumino 

Witten(WZW) class. These models are characterized by a single parameter fc, which determines the exponents. In 

this case, k turns out to be 2s for spin s. From the known behaviour of the WZW theories, Affleck showed that the 

^^ . spin spin correlations are asymptotically of the form 

un: \{s,.s,)\^\z-jr'/^'^+'\ (1) 

^D ' For any s, the field theory has a marginal operator in the renormalization group sense. This gives rise to a multiplicative 
"nI , factor la] oi [A + Bln\i — j|]^/^ in the correlation function of Eq.(^, unless the bare value of the marginal operator 
^^ I is fortuitously equal to its fbced point value [3, Q . 

"ti ■ An alternative recent line of thought has been to study explicit wave functions. In particular, the Gutzwiller wave 

Cj ' function in one dimension has been very successful. The Gutzwiller wave function takes the free Fermion determinantal 
I . wavefunction and retains all configurations with single occupancy of electrons. Thus at half filling, i.e. one electron 
''O ' per site, it yields a wave function that has only spin degrees of freedom. This insulating spin wavefunction inherits 
C , the power law correlations of the parent free Fermi wavefunction, albeit with some renormalization of the values of 
^ ' the exponents. The work of Haldane [3 and Shastry J3| and others |9| showed that the wavefunction is the exact 
• • . ground state of a long ranged s = 1/2 Heisenberg model that is in the same universality class as the Bethe chain, i.e. 
has the same correlation exponents. 

In this paper, we propose a natural extension of the Gutzwiller wavefunction to obtain wavefunctions for all s, 

Vh which (based on our numerical results) seem to have power law correlations of the form in Eq.l^. Thus we have 

constructed for wave functions the analog of the Algebraic Bethe Ansatz method for Hamiltonians: a prescription 

that yields the non-generic correlations of Eq.O automatically, i.e. without any fine tuning of parameters. This 

paper builds on ideas presented earlier in Ref. 10]. 

We start with 2s copies of the Gutzwiller wave function, and use the principle of symmetrization to produce an 
angular momentum J = s wave function, i.e. one where each lattice site has a spin s degree of freedom. Symmetrization 
is a well known procedure in angular momentum theory, where one generates the states of a spin J system by taking 
2 J copies of spin half states and projecting out all states that are not fully symmetric in the 2 J spin constituents. 
This procedure clearly generates wave functions for particles of angular momentum s. Based on the experience with 
s = 1/2, one would hope that these wavefunctions might also inherit the power law correlations of the parent free 
Fermi gas, again with some as yet undetermined renormalization of values of exponents. 



We are able to perform the symmetrization of this wave function explicitly, using the elegant formalism of spin 
coherent states. We further study the properties of this wave function using numerical techniques. When the wave 
function is squared, it can be conveniently interpreted as a statistical mechanical model of 4s parallel chains of 
logarithmically interacting particles with certain couplings, i.e. a generalization of the Wigner Dyson Coulomb 
problem. Our numerical results for s — 1,3/2,2 are presented here. Together with the analytically known results for 
s = 1/2, they support the remarkable conclusion that these wave functions provide a lattice realization of the WZW 
models with k — 2s. 

We are unable to address the issue of finding a Hamiltonian for which the states here are exact ground states, 
but have wave functions that are explicit and rather beautiful. This seems to be just the opposit situation of the 
Algebraic Bethe Ansatz models |j| for higher spin, where the Hamiltonians are relatively straightforward, but the 
wave functions are highly complicated. 

Earlier numerical work by Shastry [lOJ on the same wave function for s = 1, with relatively shorter spin chains, 
gave results that were consistent with the same exponents as found here, but lacked the resolution to determine the 
logarithmic corrections to the leading behavior. We note that spin coherent states were used for a similar mapping 
of spin wave functions to statistical mechanical problems in the work of Arovas, Auerbach and Haldane [llj as well 
as Affleck, Kennedy, Lieb and Tasaki (AKLT), who studied the s = 1 model of AKLT [l^, for which correlations 
are exponential. The calculation of the correlations is considerably easier than in the cases studied here. AKLT's 
method of solution is also a very nice application of the idea of symmetrization that we use in the present paper. 
Analytical and numerical results for the spin-half chain, with the corresponding statistical mechanical model studied 
as a function of temperature, have been presented earlier jl,l |. 

SPIN COHERENT STATES, SYMMETRIZATION AND THE GUTZWILLER WAVE FUNCTION 

Spin Coherent States 

We begin by recapitulating the salient definitions of spin coherent states. There are several related variants of these, 
and each has slightly different advantages. Let us consider a single site. The functions introduced by Radcliffe [lj| 
are defined for angular momentum J as 

\Q\ ^ g[S-cxp(#)-S+cxp(-j0)]e/2|j\ 



cos 

2J 



2./('^/2)gtan(e/2)cxp(i0)S \j\ 



= ^{'^C„}5cos2^-"(0/2)sin"(0/2)e"^|J-n). (2) 

Radcliffe also introduced a related set of states, 

\z)R^e^'^--^''^\J). (3) 

These are normalized and related to the first set by setting z -^ tan(0/2)e"^. Another set of coherent states were 
found to be very convenient |l5J for the purpose of obtaining differential operator representations of spin operators, 
e.g. S^- -^ {s — z*d/dz*) in the space of "wave functions" (^1^). These coherent states were the unnormalized states: 

\z)^e^^-\J) = Y,{^'^Cn}Hzr\J-n). (4) 

?i=0 

with |z)j7, = (1 + |zp)^-^ \z). The relationship between the wave functions il]{z*) = {z\i>), ipu{z*) =ii {z\ip) and ip{^*) 
follows from the relationship between the basis functions, wc note the relation needed later: 

V'lr^*) = cos2"'(6'/2) V(^*) with z* -^ i&n{e/2)e-"^. (5) 

The functions ip{z*) E^re polynomials in z* of maximal degree 2 J, and lend themselves to very simple "symmetrization" 
rules that we discuss in the next subsection of this paper. 

Next, we need to construct the rules for obtaining averages of variables in states. We begin by noting the diagonal 
representation of an arbitrary operator A |l6| followed by its average as: 

A= (^-^ + ^) I dna{n)\n){n\ (6) 



where the prefactor has been chosen so that o(ri) = 1 for the identity operator. To compute averages, we write 



4tt 

(2J+1) 

47r 



dn a(n)\ij(n)\^ . (?) 



Thus, given a wave function ip{^*) = (^IV-")! '^^ can find the expectation value of any operator if its a{Q) is known. 
For future reference, we note that a{il.) for the operator 5*2 is (J + l)cos6'. An important corollary that we will use 
is that if spin coherent states are constructed for every site in a lattice, the expectation value of an operator (AiBj) 
can be found by using the corresponding weight function a{fli)b{flj) in the integrals, where a{Q) and b(fl) are the 
functions for isolated sites. This will be used for the spin correlation function. 

Symmetrization 

Let us consider a simple case of two spin 2 particles with 

10) = ai/32. (8) 

Symmetrization is best understood from its action on states, so in the present case: 

S\(b)^^iaiP2+a2Pi). (9) 

On more general functions of many copies of spin half, its action is similarly defined, namely find the fully symmetric 
combination generated from a seed state and divide by the total number of generated states. We next deduce the 
rule for symmetrizing 2s copies of spin 1/2 in the space of coherent state wavef unctions. It turns out that the most 
effecient way is to work with the unnormalized coherent states EqQ) where the wave functions are just polynomials 
in z*. Let us label the 2s copies of spin 1/2 by z*{a) and the resulting spin s variables by Z, so that a coherent states 
of the direct product states and the final spin s state are: 

|{z(l), z(2) . . . z(2s)}) = n (|l/2)„ + z(a)| - l/2)„) (10) 

Q = l 

2s 

\Z) = ^f^C„}^(Z)"|2s-m). (11) 



Generic states in the direct product space and the final spin s space are represented by \4>) and [$) respectively, and 
the role of symmetrization is to map the former into the latter as 

5|0) = 10) 

10) ^ m- (12) 

The question we address is: given a state |0), what is the resulting state |$) for the spin 2s particle. The answer 
turns out to be remarkably simple in the unnormalized coherent state basis. If we denote 

/(K(l),z*(2)...z*(2s)}) = ({z*(l),z*(2)...z*(2s)}|0) 
/(K(l),z*(2)...z*(2s)}) = ({z*(l),z*(2)...z*(2s)}5|0) 

F{Z*) ^ {Z*\<^), (13) 



we find ( and demonstrate below) that 



F{Z*)^ f{Z*,Z* ...Z*). (14) 



This implies that by ignoring the distinction between the different copies of the spin 1/2 and replacing every occurrence 
of z* (a) by Z* gives us a coherent state representative of the spin s particle. This result is obvious for the special 



( symmetric) cases of / = 1 and / = ]^2;*(a), but not so obvious for other cases, since one has to rule out possible 
nontrivial dependence on the degree of the polynomial. 

We now give a brief proof of this assertion. It sufhccs to consider the general case of a polynomial of degree r, thus 



/ = z*{l)z*{2)...z*{r) 
2s- 
~2^. 



J ^ rK2srl ^ z*(^l)z*(^2)...z*(^.). (15) 



l<ll{i2...{ir<{2s) 

The state / clearly is proportional to the state deriving from {S^Y\Jz = 2s), the proportionality constant is readily 
worked out so 

m = ^^^^(^"ri2«)' (16) 

hence 

= {z*y. (17) 

The last line follows on using the commutation relations of angular momentum. 

Gutzwiller type wave functions 

The Gutzwiller wave function for a one dimensional Fermi gas at half filling is expressible in the form [g 

L 

IV'g) - ^l[^l-n,,^lHl] n iriilo) 

j=l \k\<kF 

= AA J2 e^"^.'-^\l[sm'inir,-rk)/L)]S-S-^...S-^ \L/2) (18) 

l<ri{. ..{r„<L k{j 

with n = L/2. Here cj. is a creation operator with wavevector k in the original Fermionic representation, and the 
[1 — nj^-^rij^i] factors ensure no double occupancy for any site j. S~ are the spin lowering operators in the spin 
representation that is equivalent at half filling. Af is the normalization that we will not specify till the end, since it 
cancels out in the evaluation of the correlations. For simplicity we have confined our considerations to the case of L/2 
overturned spins, so we are dealing with a global singlet wave function, made up of L spin 1/2 particles. 

We next consider 2s copies of this wave function, and project into the spin s sector at each site. In view of the 
discussion in the last section, this is most easily done with the coherent state notation, so the product wave function 
is written down directly as 

^{Zl,Z;,...Zl) ^ M[ J2 e^''^S.'-^>[nsin'('r(r,-r,)/L)]Z;z;...Z;^^j2^ (19) 



I L/2' 

l<ri{...{r_L<L k<j 



(20) 



To obtain this result, we wrote \^g) in terms of z*{a), multiplied 2s copies of this, and then symmetrized the 
wavefunction by dropping the distinction between the different copies or replicas. To reconstruct the wave function 
in the angular basis we use Eq.lj^l and write 

L 

^{ni,n;...ni) = l[cos^%e,/2)^{z:^,z:^...z:j. (21) 

In order to make this more tractable, we introduce "occupation numbers" p" which determine whether we get a 
cos{9j/2) or sin(6'j72) factor at a given site j in a particular " replica" a. Thus at each site we get a factor of 

2s 

factor = Yl [cos(0j/2)(l - p'^) + sin(0j/2)e-*+*''V5'] 



2s 



= cos^''{ej/2)exp[^ppn{tan{ej/2)+iTrj~i(l)j}]. (22) 



Q = l 

We thus write the wave function as 



i=i {p°=o,i} 






(23) 
The sum over the occupancy integers p" is constrained to obey ^„ p" — L/2 for each repHca a. In the next section 
we continue the discussion of expectation values and correlation functions, which require taking the modulus square 
of this wave function. 

COUPLED ISING REPRESENTATION 

The probability density is |^({il| . . . il^Dl^i so ^'^ need to multiply ^ in Ea. (|23|l with its complex conjugate, 
leading to 4s rephcas in all. The correlation function (s|s|) can be found by calculating (s + l)'^{cos9j cos9k) with 
this probability density, as noted after Eq(|7|). From the fact that |^) is a singlet state, this correlation function suffices 
to determine all components (s'^s^). The angular variables 9 and </> can be integrated over, leaving only the pj's. The 
problem reduces to one of interacting lattice gas particles. 

Integration over the azimuthal angle (j)j can be done at each site j, and gives a constraint that ^^ p" — ^^, p" = 
s — TOj, where a and a' refer to replicas in Vt and \1/* respectively. Note that the oscillating phase factor can be 
dropped in view of the constraint from the azimuthal integration. 

The integral over 9j is next performed: 

^(sin|)2-^"(cos 1)^^+2" sin e,d0, = 2^-^±0^^^ = 2W.im), (24) 

with the function Ws{m) is defined by this equation. For the spin autocorrelation function {s^s^}, we have to insert 
an extra factor of (s + 1) cos9 in the 9 integrals at j and k. It is easy to verify that 

(s + 1) / (sin ^)2^-2"(cos ^)2«+2™ cos9j sm9jd9j = 2mWs{m). (25) 

Jo 2 2 

The final result is like a classical partition function: 
Z^\^i{ni...ni})\'=N- E expEp>fclnsin2(7r(j-fc)/L) + Ep">fc'lnsin2(^(j-fc)/L)]nW^.(™.) (26) 

where a and a' refer to replicas in \I/ and \I/* respectively, rrij — s — J^a P? ^-'^'^ ^^'^ ^^^^ ^^ ^^"^ partition function Z 
is subject to the constraints 

T.P? = E^"'-V2 

3 3 

E^i - Epf- (27) 

a a' 

When calculating the spin autocorrelation function (SjSl), the summation in Ea. (|26|l for Z is evaluated with an extra 
factor of rrijiTik. 

It is convenient to change variables to a" = 2p" — 1 and cr" = 1 — 2p" . We then have 4s coupled Ising chains (2s 
from ^ and 2s from '5*). The constraint Ea. (|27^ is equivalent to the condition that the sum of the Ising spins (not 
to be confused with the original quantum spins) across the different chains at any site must be zero. In addition, we 
impose the condition that the sum of the spins along any chain must be zero. With this condition, the interactions 
along any chain are 

E ^^^^ lnsin2(vr(j - k)/L) = const + E J<^." Insin2(7r(j - k)/L). (28) 

j<k j<k 



The interactions along any chain are antiferromagnetic and logarithmic. The interaction across the different chains 
at any site occurs, apart from the constraint, through the Wsijn) factor. 

For the case of s = 1, there is an alternative form of the partition function that is more convenient. With the 
constraint that the spins on the four chains at any site must add up to zero, there are six possible configurations 
at any site: (1, 1, —1, —1), (1, —1, 1, —1), (1, —1, —1, 1), and the mirror images of these three. The magnetization m 
is 1 for the first and zero for the next two configurations. A remarkable simplification occurs once we note that 
^^ (t"(t^F(j — k) is equal to zero for any function F unless the configurations at the sites j and k are either identical 
or mirror images. (In the sum over a, the two replicas from ips and the two from ip* are included.) Thus instead of 
four coupled Ising chains, the problem reduces to a six state Potts model on one chain. The states are labelled by 
q = ±1, 2, 3. The configuration (1, 1, 1, 1) is labelled with q = 1, so that rrij = ±1 when qj = ±1, and rrij — when 
qj — ±2, ±3. The partition function is 

Z ex ^ exp[^(<5,^,,, - <5,^,_,J lnsin2(^(j - k)/L)] J]!! + 5,|,i)- (29) 

qi---qL j<k j 

The last factor comes from the fact that Ws(±l) = 2Ws{0). Two sites only interact with each other if their g's are 
identical or opposite. The condition that the total magnetization for each of the four original chains must be zero 
reduces to the statement that 

for any n. 

Potts models can similarly be constructed for s > 1, although they are more complicated: two sites j and k 
interact even when qj ^ iqk, and the condition from the total magnetization is weaker than Eq.ljSOJ. The numerical 
simulations reported in the next section were conducted with both the coupled Ising and the Potts representations. 

NUMERICAL RESULTS 

Monte Carlo simulations were performed on the Ising and Potts representations given by Eas. (|26|l and (|29|l . For 
the Ising case, the constraint that the total spin along any chain and for all replicas at a site must be zero prevents 
using single spin flip dynamics. In a Monte Carlo move, two replicas (a, /3) on two adjacent sites (j, k) were chosen at 
random. If the Ising spins at these four locations satisfy a'^ — a^, a^ — cr^ and a'^ y^ cr^, it is possible to flip all the 
four spins simultaneously without violating the constraint. The ratio of the probability of the flipped configuration to 
the probability of the unflipped configuration is calculated, with appropriate factors of Ws (m) included if the move 
would change the magnetization at the sites, i.e. if a < 2s < /3 or /3 < 2s < a. The move is accepted or rejected using 
the standard Metropolis criterion. 

For the Potts representation, two adjacent sites (j, fc) were chosen at random, and {qj,qk) were attempted to be 
changed. For the case of s = 1, Ea. (|Sn)l required that if qj + qk ^ 0, the only possible move was to exchange them, 
while if qj + q^ = 0, one could attempt to replace them with either of the other two pairs of g- values. For s > 1, a 
table was constructed at the beginning of the numerical simulation. For any pair of q- values (gi, (72), the table listed 
all pairs {q'l, q'2) that {qi, 92) could change to, while respecting the magnetization constraints on the underlying Ising 
spins. In any Monte Carlo step, this table was used to randomly select an allowable move to attempt. For both 
the Ising and the Potts representations, the long range logarithmic interaction down the chains made calculating the 
probability of an attempted move an 0{L) long calculation for a chain of length L. As a result, very large values of 
L could not be simulated. 

In both cases, error bars on the measured correlation function were estimated by taking blocks of 30000 readings, 
calculating the average correlation function within each block and then the inter-block variance. Even if the individual 
readings are taken too frequently and are therefore correlated, this procedure should be reliable so long as the blocks 
are sufficiently large to be uncorrelated. It is also useful to compare the variance of the block averages to the variance of 
the individual readings within a block. The latter would be 30000 times the former if the readings were uncorrelated. 
For our simulations, the actual ratio ranged from about 30000 to 300, confirming that in all cases the blocks are 
uncorrelated even when the individual readings are not. 

Figure^shows the spin autocorrelation function C{r) — {—ly (sjSj+r) as a function of r for different L, for s = 1. 
The antiferromagnetic interaction down the chains causes the oscillatory (—1)*^ factor. It is clear from the figure that 
the correlation function does not decay exponentially. As seen in the figure, the autocorrelation function seems to 
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FIG. 1: Log-log plot of C{r), the magnitude of the spin autocorrelation function, as a function of r for s — 1. System sizes 
ranging from L = 128 to L = 2048 are plotted. The error bars are comparable to the point size. The data is consistent with a 
power law form, with exponent —0.68. Since logarithmic corrections are expected, a scaling collapse is not shown. 
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FIG. 2: Plot of C^{r)r^'^ as a function of Inr for s = 1. This is predicted to be a straight line, with finite size corrections. 
System sizes from L = 128 to L = 2048 are plotted. The error bars are comparable to the point size. 



decay as l/r^ with /i : 
functional form 



0.68. However, motivated by the analytical considerations discussed in Section 1, we try the 



C(r) 



[A + Blnr]^^^ — 



(31) 



with /i = 0.75. This is because the power-law part of C{r) should scale as ^ l/r'' with ^ = 3/(2s-|-2) — 0.75 for s = 1, 
with possible logarithmic corrections from the marginally irrelevant operator. Accordingly, Figure |21 plots C^(r)r^'^, 
which should be a linear function of Inr. This expectation is borne out by the plot. Logarithmic dependences are 
known to be hard to distinguish from weak power laws, and the plots cannot be used to choose between Ea. (|31|l and 
the pure power law decay of Figure ^ (It would be even harder for the data to discriminate between more subtle 
differences, e.g. Ea. (|31|l with different exponents to the logarithmic term in the numerator.) However, from FigureEl 
the data is certainly in agreement with the analytical expectation. 

Figure |2| and Figure 0] are the counterparts of Figure Q] and Figure [21 respectively for s = 3/2. Based on the 
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FIG. 3: Log-log plot of C(r), the magnitude of the spin autocorrelation function, as a function of r for s = 3/2. System sizes 
ranging from L — 128 to L = 2048 are plotted. The error bars are comparable to the point size. The data is consistent with a 
power law form, with exponent —0.51. 
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FIG. 4: Plot of C'^{r)r^'^ as a function of Inr for s — 3/2. System sizes from L — 128 to L = 2048 are plotted. The error bars 
are comparable to the point size. 



analytical prediction, Figure0]plots C'^{r)r^''^ as a function of Inr, since the power law decay part of C{r) should have 
an exponent oi ^ — 3/(2s + 2) — 0.6. Likewise, FigureEland FigureElare for s = 2, with C^{r)r^-^ plotted in FigureEl 
Although the results for s = 3/2 and s = 2 are not as clear as those for s = 1, we see that C{r) definitely does not 
decay exponentially for s — 2, or as ~ 1/r for s = 3/2, the generic behavior expected for integer and half-integer spin 
chains respectively. The decay of C{r) is consistent with the expectation that each s corresponds to a WZW model 
of order k = 2s. The effective /x's from Figures 01 and El are 0.51 and 0.40 respectively. 

Since there is no clear linear region in Figure^l we also plot the data in a manner that eliminates finite size effects. 
For finite system size L, the leading effect on C{r) is to change its form from Ea. H31(l to 



C(r) 



[A + Blnr]^^^^F{r/L) 



(32) 



where F is an unknown function. Since ^ — 0.5 for s = 2, we plot C^{r)r for fixed r/L as a function of InL. The 
result should be a straight line, with slope BF'^{r/L). The results are shown in Figure [3 For all the values of r/L 
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FIG. 5: Log-log plot of C{r), the magnitude of the spin autocorrelation function, as a function of r for s — 2. System sizes 
ranging from L — 64 to L — 1024 are plotted. The error bars are comparable to the point size. The data is consistent with a 
power law form, with exponent —0.40. 
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FIG. 6: Plot of C'^{r)r^''^ as a function of Inr for s — 2. System sizes from L = 64 to L = 1024 are plotted. The error bars are 
comparable to the point size. 

shown, there is an upward curvature to the plots. This is presumably due to subleading corrections to the scaling 
form, since the leading correction is large enough to change the effective /i in Ea. (|31|l from 0.5 to 0.4. 



CONCLUSION 



In this paper, we have proposed a generalization of the Gutzwiller wavefunction for spin i chains that yields simple 
wavefunctions for spin chains with s > ^ Remarkably, the spin spin correlation functions for these wavefunctions 
have the same power law decay as Wess Zumino Witten models with k = 2s, in contrast to the generic expectation of 
k = \ for all half-integer spin and an exponential decay for integer spin. This result was obtained through numerical 
simulations for s = 1, | and 2, after mapping the model to a classical statistical mechanical model with long range 
interacting spin chains. 

In summary, we have taken free fermionic wave functions and found a way of projecting them in a fashion that 
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FIG. 7: Plot of C {r)r for fixed r/L as a function of logj L, for s — 2. Four different values of r/L are shown. The vertical 
bars are the error bars. There is a noticeable odd-even effect for r — L/8 and r — 3L/8 for L = 8. For all the values of r/L, 
the plots curve slightly upward. 

yields the WZW theory exponents, a possibility that has been presaged in Ref(|l3)- It remains to be seen if there is 
a systematic way of finding Hamiltonians for which the wave functions presented here are ground states. 
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